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ABSTRACT 

We present transit photometry of three exoplanets, TrES-4b, HAT-P-3b, and WASP-12b, allowing for refined 
estimates of the systems' parameters. TrES-4b and WASP-12b were confirmed to be "bloated" planets, with 
radii of L706± 0.056 Rjup and L736± 0.092 Rjup, respectively. These planets are too large to be explained with 
standard models of gas giant planets. In contrast, HAT-P-3b has a radius of 0.827 ±0.055 Rjup, smaller than 
a pure hydrogen-helium planet and indicative of a highly metal-enriched composition. Analyses of the transit 
timings revealed no significant departures from strict periodicity. For TrES-4, our relatively recent observations 
allow for improvement in the orbital ephemerides, which is useful for planning future observations. 
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1. INTRODUCTION 

A puzzling feature of the hot Jupiters is that many of them 
have radii that are either larger or smaller than one would 
have guessed prior to the discovery of this class of objects. 
This "radius anomaly problem" has been present sin ce the 
first transiting^ planet was discovered by OiarbonneaiTetari 
d^OO) and iHenry et al.l (l200(l. and stiU has no universally 
acknowledged resolution. lFortney & Nettelmannl(12010 ) have 
reviewed many of the proposed solutions, and even in the 
short time since thei r review several ot her theorie s have been 

fjposed (see, e.g., I Perna et 311120101: ^atvgin & StevensonI 
10). The small size of some planets can be explained as a 
consequence of heavy-metal enrichment, beyond the enrich- 
ment factors of Jupiter and Saturn and comparable to those of 
Uranus and Neptune. As for the larger radii, possible expla- 
nations include tidal friction, unexpected atmospheric prop- 
erties, and resistive heating from electrical currents driven by 
star-planet interactions. 

This paper presents follow-up observations of three exo- 
planets that were found to have anomalous radii. As in pre- 
vious papers in this series, the purpose of the observations 
was to refine the system parameters (thereby checking on the 
magnitude of the radius problem) and to check for any transit 
timing anomalies that might be caused by additional gravitat- 
ing bodies in the system. 

Two of our targets are among t he most "bloated" planet s 
known. TrES-4 was discovered bv iMandushev et akl (120071) . 
and the two high-precision light curve s that accompanied 
the discove ry paper were reanalyzed by [Torres et"ari (t2008l l 
and ISozze tti et al. (2009). Here, we present 5 new light 
curves for this system. We also present 2 new l ight curves 
for WA SP- 12, a planet that was discovered by iHebb et al.l 
( 120091) and for which occultation photometry has been used to 
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characterize the planet's atmosp here and orbit (ICampo et al.l 
'2010'; Madhusudhan elaLmoTl). Our third target, HAT-P-3 
(Torres et al. 2007), i s in the o pposite category of planets that 
are "too small." Gib son et al.l (12010) have published 7 high- 
quality light curves of the system. We present six new light 
curves, and provide independent estimates of the planetary 
and stellar parameters. 

2. OBSERVATIONS AND DATA REDUCTION 

Almost all the observations were conducted at the Fred 
Lawrence Whipple Observatory (FLWO) located on Mt. Hop- 
kins, Arizona, using the 1 .2m telescope and KeplerCam de- 
tector The KeplerCam is a 4096^ CCD with a field of view 
of 23. '1 X 23.'1. The pixels were binned 2 x 2 on the chip for 
faster readout. The binned pixels subtend 0."68 on a side. Ob- 
servations were made through Sloan / and z filters. One of 
the WASP- 12 transits was observed with the Nordic Optical 
Telescope (NOT) located in the Canary Islands, using the AL- 
FOSC detector The ALFOSC detector is a 2048 ^ CCD with 
a field of view of 6. '4 x 6. '4, corresponding to 0"19 per pixel. 
The observation was made through a Johnson V filter On 
each night we attempted to observe the entire transit, with at 
least an hour before ingress and an hour after egress, but the 
weather did not always cooperate. 

We performed overscan correction, trimming, bias subtrac- 
tion and flat-field division with IRAF0. To generate the light 
curves, we performed aperture photometry on the target star 
and all the comparison stars with similar brightnesses to the 
target star (within about a factor of two). We tried many dif- 
ferent choices for the photometric aperture and found, unsur- 
prisingly, that the best aperture diameter was approximately 
twice the full width at half maximum (FWHM) of the star 
A comparison signal was formed from the weighted average 
of the flux histories of the comparison stars. The weights 
were chosen to minimize the out-of-transit (OOT) noise level. 
Some comparison stars that did not seem to provide a good 
correction were rejected. In general, 6-8 comparison stars 
were used to generate the final comparison signal. The tar- 
get star's flux history was divided by the comparison signal 

^ IRAF is distributed by the National Optical Astronomy Observatory, 
which is operated by the Association of Universities for Research in Astron- 
omy (AURA) under cooperative agreement with the National Science Foun- 
dation. 
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Fig. 1 .— Relative photometry of TrES-4 in the /-band. From top to bottom, the observing dates are 2008 Jun 10, 2009 Apr 1, 2009 May 3, 2010 Apr 27 and 
2010 May 4. See Table 1 for the cadence and rms residual of each light curve. The bottom plot is a composite light curve averaged into 3 min bins. 



and then multiplied by a constant to give a mean flux of unity 
outside of the transit. 

Table 1 is a journal of all our observations, including those 
that were spoiled by bad weather. The light curves are dis- 
played in Figures 1-4, after having been corrected for differ- 
ential extinction. Section 3 explains how this correction was 
applied. The airmass-corrected data are given in electronic 
form in Tables 2-4. 



3. DETERMINATION OF SYSTEM PARAMETERS 

Our techniques for light-curve modeling and parameter 
estimation are similar to those employed in previous pa- 
pers in this series (see, e.g., Holman et al. 2006; Winn et al. 
I2007ah. The basis f or the light-curve model was the formula 
of iMandel & Ago! (l200l . assuming quadratic limb darken- 
ing and a circular orbit. The set of model parameters in- 
cluded the planet-to-star radius ratio (Rp/Ri,), the stellar ra- 
dius in units of orbital distance (R^/a), the impact parameter 
(b = acosi/R^, where / is the orbital inclination), the time 
of conjunction for each individual transit (Tc), and the limb- 
darkening parameters U[ and U2- We also fitted for two param- 
eters ( Amo, k.) specifying a correction for differential extinc- 
tion, 

Amcor = Amobs + Amo + k^Z, (1) 



where z is the airmass, Amobs is the observed magnitude, and 
Amcor is the corrected magnitude that to be compared to the 
idealized transit model. 

Since the data are not precise enough to determine both of 
the limb-darkening parameters, we followed the suggestion of 
iPall (I2OO81) to form uncorrected linear combinations of those 
parameters. We allowed the well-constrained combination to 
be a free parameter and held the poorly-constrained combi- 
nation fixed at a tabulated value for a star of the appropriate 
type. In Pal's notation, the rotation angle was taken to be 
near 37° in all cases. To determine the tabulated values we 
used a program kindly pro vided by J. S outhworth to query 
and interpolate the tables of IClareli (120041) 11 

For parameter estimation, we used a Markov Chain Monte 
Carlo (MCMC) algorithm to sample from the posterior prob- 
ability distribution, employing the Metropolis-Hastings jump 
criterion and Gibbs sampling. Uniform priors were adopted 
for all parameters, and the likelihood was taken to be 
exp(-x^/2) with 




' The interpolated values (ui.m) for the cases of TrES-4 i, HAT-P-3 
i, HAT-P-3 z, WASP-12 z, and WASP-12 V are (0.20,0.37), (0.38,0.27), 
(0.30,0.29), (0.14,0.36), (0.36,0.35), respectively 
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Fig. 2. — Relative photometry of HAT-P-3 in the /-band. From top to bottom, the observing dates are 2008 Mar 8, 2008 Apr 6 and 2008 May 5. See Table 1 for 
the cadence and rms residual of each light curve. The bottom plot is a composite light curve averaged into 3 min bins. 



where fobs.ij is the jth data point from the ith Ught curve, 
/caici.j is the calculated light curve based on the current pa- 
rameters, and aij is the uncertainty associated with fobs.ij- 
All the light curves for a given planet were fitted simultane- 
ously. The uncertainties were determined in a two-step pro- 
cess. First, the standard deviation ct, of the OOT data was 
determined for each light curve. In a few cases, the pre- 
transit and post-transit noise levels were very different (due 
to different airmasses); in these cases the starting point was 
a function a, j that interpolated linearly between the two dif- 
fering noise levels. Second, the preceding uncertainty esti- 
mates were multiplied by a correction factor (3 > I intended 
to account for time-correlated noise. The /3 factor was deter- 
mined with the "time-averaging" proce dure (iPont et aTll2006t : 
iWinn et al.l 2008': 'Ca rter & Winn|2009l) . using bin sizes brack- 
eting the ingress/egress duration by a factor of 2. The values 
of /3 are given in Table [T] 

The starting point for each Markov chain was determined 
by minimizing x^, and then perturbing those parameters by 
Gaussian random numbers with a standard deviation of 10a, 
where a is the rough uncertainty estimate returned by the 
least-squares fit. We ran several test chains to establish the ap- 
propriate jump sizes, giving acceptance rates near 40%. Then 
we ran 4-5 chains each with 10* links, ignored the initial 
20% of each chain, and ensured convergence according to the 



Gelman-Rubin statistic dOelman & Rubinlll992h . The quoted 
value for each parameter is the median of the one-dimensional 
marginalized posterior, and the quoted uncertainty interval en- 
closes 68.3% of the probability (ranging from the 15.85% to 
the 84.15% levels of the cumulative probability distribution). 

4. RESULTS 

The results for the model parameters, and various derived 
parameters of interest, are given in Tables|5]|6]and|7] The next 
subsection explains how the stellar and plantetary dimensions 
were calculated from the combination of light-curve param- 
eters and stellar-evolutionary models. This is followed by a 
subsection presenting an examination of the transit times. 

4.1. The stellar and planetary radii 

The stellar and planetary radii and masses cannot be deter- 
mined from transit parameters alone. The route we followed 
to determining these dimensions was to set the mass scale by 
using an estimated stellar mass and the observed semi- 
amplitude A!* of the star's radial-velocity orbit. The stellar 
mass is itself estimated by using stellar-evolutionary models 
with inputs from the observed spectral parameters, as well 
as the mean density that is calculated from the light-curve 
parameters. For the re levant form ulas and discussion see 
ISozzetti et all ( l2007h or IWinnI (|20!0). We used previously 
measured values of Ki,, documented in Tables |5] |6] and Q 
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Fig. 3. — Relative photometry of HAT-P-3 in the z-band. From top to bottom, the observing dates are 2009 Mar 14, 2009 Mai- 20 and 2009 Apr 15. See Table 
1 for the cadence and rms residual of each light curve. The bottom plot is a composite light curve averaged into 3 min bins. 



For the evolutionary models, we used the Yonsei-Yale (Y^) 
isochrones (Yi et al. 2001). 

The Y^ isochrones can be thought of as an algorithm that 
takes as input the age, metallicity, mass, and concentration 
of a-elements of a star, and returns the star's temperature, 
mass, density, and other properties. We interpolated the Y^ 
isochrones in age from 0.1 to 14 Gyr in steps of 0.1 Gyr, 
and in metallicity from -0.20 to 0.58 dex in steps of 0.02 
dex. Then we used linear interpolation to create a 4 x finer 
mass sampling for each metallicity and age. We assumed 
the concentration of a-elements to be solar To each model 
star in the resulting Y^ grid, we assigned a likelihood based 
on the measured metallicity Z and effective temperature T^ff 
(taken from the literature) as well as the stellar mean density 
/3,^ determin e d sole ly from the transit parameters. Following 
iCarter et al] (12009b . the likelihood was taken to be propor- 
tional to «exp(-x^/2), where 

^ / Z-Z^\ ' ^ / reft- ref,ob. y ^ j^^y ^ 
\ (Tz J \ (J%„ J VP*.obs/ 

and n is the number density of stars as a function of mass, 
according to Salpeter's law with exponent -1.35. The effect 
of multiplying by n is to set a prior so that extremely rare and 
short-lived stars are disfavored, even though they might pro- 
vide a good fit. (The effect of this prior was generally small.) 



Finally, the "best-fitting" values and uncertainties were com- 
puted from the appropriate likelihood-weighted integrals in 
the space of model stars. 

With the stellar mass thereby determined, it is possible to 
compute the other dimensions Rp,Ri,, and Mp. The results 
are given in Tables |5] |6] and [T] We note that this proce- 
dure does not take into account any uncertainty in the Y^ 
isochrones themselves and, therefore, is subject to system- 
atic errors that pr obably amount to a few percent (see, e.g., 
iTorres et al.ll200l) . 

4.2. Transit times and revised ephemerides 

We analyzed the new transit times in conjunction with pre- 
viously published midtransit times, to seek evidence for sig- 
nificant discrepancies from strict periodicity and refine the 
ephemerides to allow for accurate prediction of future events. 
The new transit times are given in Tables |8] |9] and [TO] with 
uncertainties determined by the MCMC analysis described in 
Section 3. Previously published midtransit times were taken 
from Gibson et al. (2010), Hebb et al. (2009), Sozzet ti et all 
(f2009.) . . Torres et al.i (i2007i) . and iMandushev et al.. (,2007Ik 

The transit times for each system were fitted with a linear 
function, 

Tc = To + EP (4) 
where E is an integer (the epoch), P is the period, and Tq is a 
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Fig. 4. — Relative photometry of WASP-12 in the z- and V-bands. The top light curve is based on z-band observations on 2009 Jan 08. See Table 1 for the 
cadence and rms residual of each light curve. The bottom light curve is based on V-band observations on 2009 Dec 06. 



particular reference time. The best-fitting values of P and To 
were determined by linear regression. Figures 5-7 show the 
residuals to the fits, and the captions specify the minimum 
and the number of degrees of freedom. In no case is there 
clear evidence of timing anomalies. For TrES-4, there is for- 
mally only a 6% chance of obtaining such a large with 
random Gaussian errors, but we do not deem this significant 
enough to warrant special attention. 

The c ase of WASP -12 required somewhat special treatment 
because iHebb et all (12009) did not report individual midtran- 
sit times, but rather a consensus reference time based on ob- 
servations of multiple transits. We included their qu oted refer- 
ence time as a single data point. iCampo et al.l (12010 ) provided 
many transit times obtained over several years, but most of the 
data were from amateur observers and were not presented in 
detail or evaluated critically. For this reasons we did not in- 
clude them; however, as a consequence, there are not many 
points in our fit. 

The results for P and Tq are given in Tables 5-7. They are 
based on a fit in which all of the uncertainties of the transit 
times were rescaled by a common factor to give x^/^dof = 
1 . Our intention is to provide conservative error estimates to 
allow for planning of future observations. The uncertainties 
on the individual transit times given in Tables 8-10 were not 
rescaled in this way, nor were the error bars that are plotted in 



Figures 5-7. 
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Fig. 5. — Timing residuals for TrES-4. The data presented in this paper are 
labeled with solid circles (complete transits) and solid trianges (partial tran- 
sits). Open circles represent transits observed by Mandushev et al., t2007i) . 
The best fit gives - 12.1 with 6 degrees of freedom. The probability of 
obtaining a higher v/ith random Gaussian data points is about 6%. 



5. SUMMARY AND DISCUSSION 

We have presented new photometry and new analyses of the 
transiting exoplanets TrES-4b, HAT-P-3b, and WASP-12b. 
Whereas the discovery papers reporting TrES-4b and HAT-P- 
3b included only a few high-precision light curves, our analy- 
ses are based on 5-6 such datasets. Likewise, the WASP-12b 
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Fig. 6. — Timing residuals for HAT-P-3. The data presented in tliis paper 
are labeled with solid circles (complete transits) and solid trianges (partial 
transits). Open circles represent transits observed by Gibson et al. 12010); 
[Torres et al.l <2007l) . The best fit gives = 10.4 with 12 degrees of freedom. 
The probability of obtaining a higher with random Gaussian data points is 
about 60%. 




Epoch 

Fig. 7. — Timing residuals for WASP-12. Solid circles are complete tran- 
sits that we observed, and the open circle is the reference transit time (derived 
from observations of multiple events) quoted bv .Hebb et al. (2009.) . The best 
fit gives = 2.0 with 1 degree of freedom. 



discovery paper featured only one high-precision light curve, 
to which we have added two. We have applied consistent and 
conservative procedures for parameter estimation, including 
an accounting for uncertainties in the limb-darkening law and 
due to time-correlated noise, as well as linkage between the 
light curve parameters and stellar-evolutionary models, that 
were not always applied by previous authors. In the past 
these efforts have occasionally led to significant revisions of 
the planetary dimensions (see, e.g. Winn et al. 2007b, 2008). 

In the present case our results are in agreement with the pre- 
viously reported results. All of the TrES-4 parameters agree 



to within 2<7 with the results reported bv lTorres et al.l (12008b . 
Two of the most important parameters, the mass and radius of 
the planet, agree to within la. Our HAT- P-3 param eters agree 
to within 2a with those reported by Gibs on et al.l (2010), and 
the planetary mass and radius agree to within 1. 3a. Our 
W ASP-12 parameters agree to within 2a with those reported 
bv lHebbet all (l2009l) . and the mass and radius of the planet 
agree to within la. 

Our transit ephemeri s for TrES-4 agrees with that of 
iMandushev et alj ( 120071) . and the refined orbital period is 
about 20 times more precise. Our eph emeris for HAT-P- 

3 agrees with that of iGibson et al.l ("2010), and is 2-3 times 
m ore precise. Our e phemeris for WASP-12 agrees with that 
of iHebb et alj (120091) . and is of comparable precision, despite 
being based on only a few well-documented data points. 

The reason that TrES-4, HAT-P-3, and WASP-12 ai-e of 
particular interest is because their measured dimensions do 
not agree with standard models of gas giant planets. TrES- 

4 and WASP-12 are heavily bloated, with radii too large for 
their masses , while HAT-P-3 is to o small. With reference to 
the tables of iFortnev et"aLl (l2007h . TrES-4 and WASP-12 are 
incompatible with pure hydrogen-helium giant planets at the 
lOcr and 5a levels, respectively. Enhancing these planets with 
metals to the degree of Jupiter or Saturn would only make the 
problem worse. HAT-P-3, in contrast, is compatible with the 
tabulated models if it is endowed with approximately lOOM^ 
of heavy elements. Our results do not change these interpreta- 
tions of the three systems we have studied. Rather, our results 
lend more confidence to the claims that the dimensions of the 
planets are anomalous, and merit attention by theoreticians 
who seek to solve the radius anomaly problem. 
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TABLE 1 
Journal of Observations 



Date 


Target 


Filter 


Cadence 


RMS 


Red noise factor 


Notes 


[UT] 






r [min"'] 


fj 


/3 




08 Mar 2008 


HAT-P 3 


i 


1.54 


0.0018 


1.39 




06 Apr 2008 


HAT-P-3 


i 


2.50 


0.0019 


1.48 




05 May 2008 


HAT-P-3 


i 


2.50 


0.0036 


1.00 


Incomplete phase coverage 


18 Jan 2009 


HAT-P-3 


z 


0.58 


0.55 


_ 


Bad weather" 


14 Mar 2009 


HAT-P-3 


z 


0.81 


0.0010 


1.31 




20 Mar 2009 


HAT-P-3 


z 


0.67 


0.0014 


1.24 


Incomplete phase coverage 


15 Apr 2009 


HAT-P-3 


z 


0.81 


0.0012 


1.05 




18 Apr 2009 


HAT-P-3 


z 


0.81 


0.0013 


- 


Bad weather'* 


30 Jan 2010 


HAT-P-3 


z 


0.58 


0.0029 


- 


Very incomplete phase coverage'* 


10 Jun 2008 


TrES-4 


i 


0.81 


0.0013 


1.06 


Incomplete phase coverage 


01 Apr 2009 


TrES-4 


i 


0.81 


0.0012 


1.00 




03 May 2009 


TrES-4 


i 


0.67 


0.0009 


1.42 


Incomplete phase coverage 


28 May 2009 


TrES-4 


i 


0.71 


0.0011 


_ 


Very incomplete phase coverage" 


06 Jul 2009 


TrES-4 


i 


0.67 


0.0013 


_ 


Very incomplete phase coverage" 


27 Apr 20 10 


TrES-4 


i 


0.45 


0.0016 


1.00 




04 May 2010 


TrES-4 


i 


0.67 


0.0012 


1.15 


Partial transit only 


16 Jun 2010 


TrES-4 




0.81 


0.0011 




Very incomplete phase coverage" 


18 Dec 2008 


WASP- 12 




0.80 


0.0020 




Bad weather" 


08 Jan 2009 


WASP- 12 




1.50 


0.0018 


1.48 




18 Jan 2009 


WASP- 12 




0.67 


0.0019 




Bad weather" 


19 Jan 2009 


WASP- 12 


z 


1.00 


0.0011 




Bad weather" 


07 Mar 2009 


WASP-12 


z 


0.80 


0.0012 




Bad weather^ 


06 Dec 2009 


WASP-12 


V 


6.33 


0.0011 


1.57 


Observed using the NOT 


12 Jan 2010 


WASP-12 


z 


1.00 


0.0020 




Bad weather" 


24 Jan 2010 


WASP-12 


z 


0.50 


0.0031 




Bad weather" 


25 Jan 2010 


WASP-12 


z 


0.50 


0.0023 




Bad weather" 


18 Feb 2010 


WASP-12 


z 


0.40 


0.0018 




Bad weather" 


01 Mar 2010 


WASP-12 


z 


0.80 


0.0049 




Bad weather" 



" Data not used in calculations. 



TABLE 2 
Photometry of TrES-4 (Excerpt) 



Filter HJDurc Relative flux 


Airmass 


i 2454628.84123771 1.00133657 
i 2454628.84212893 0.99907639 
i 2454628.84384188 0.99966428 


1.0055 
1.0058 

1.0062 


Note. — The time-stamp represents the UT-based Helio- 
centric Julian Date at midexposure. We intend for the rest of 
this table to be available online. 


TABLE 3 

PllOTOMinRY ()!■ H.\T-P-3 (EXC'LiRP 




Filter HJDutc Relative flux 


Airmass 


i 2454534.74665524 1.00252886 
i 2454534.74709502 0.99836326 
i 2454534.74755801 1.00316729 


1.5209 
1.5175 
1.5142 



Note. — The time-stamp represents the UT-based Helio- 
centric Julian Date at midexposure. We intend for the rest of 
this table to be available online. 



TABLE 4 

Photometry of WASP-12 (Excerpt) 



Filter 


HJDutc 


Relative flux 


Airmass 


z 


2454840.62025635 


1.00042427 


1.5419 


z 


2454840.62085820 


1.00001333 


1.5362 


z 


2454840.62148318 


1.00215610 


1.5304 



NOTli. — The time-stamp represents the UT-based HeUo- 
centric Juhan Date at midexposure. We intend for the rest of 
this table to be available online. 
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TABLE 5 
System Parameters for TrES-4 



Parameter 


Symbol 


Value 


68.3% Conf. Limits 


Comment" 


Transit parameters: 










Orbital period [d] 


P 


3.5539268 


±0.0000032 


LC 


Midtransit time [HJDujc] 




z454z3U. 90524 


±0.00062 


LC 


Planet-to-star radius ratio 


J? / J? 

I 




±u.uuu/o 


T 


Orbital inclination [deg] 


oZ.ol 


±\}.i 1 


T 


Scaled semimajor axis 


a/R-, 


6.08 


±0.16 


LC 


Transit impact parameter 


b 


0.761 


±0.018 


LC 


Transit duration [hr] 




3.567 


±0.037 


LC 


Transit ingress or egress duration [hr] 




0.692 


±0.047 


LC 


Other orbital parameters: 










Orbital eccentricity 


e 







Adopted 


Velocity semiamplitude [m s"' ] 


K 


97.4 


±7.2 


K 


Planet-to-star mass ratio 


Mp/M^ 


0.000631 


±0.000047 


LC + M + T+Y^ + K 


Semimajor axis [AU] 


a 


0.05084 


±0.00050 


LC + M + T + Y2 


Stellar parameters: 










Mass [Mq] 




1.388 


±0.042 


LC + M + T + Y2 


Radius [Rq] 




1.798 


±0.052 


LC + M + T + Y2 


Mean density [g cm"^] 


P* 


0.337 


±0.027 


LC 


Errective temperature [K] 




6200 


±75 


T 


Projected rotation rate [km s"'] 


vsiniV 


9.5 


±1.0 


PRR 


Surface gravity [cgs] 


lege* 


4.071 


±0.024 


LC + M + T + Y^ 


Metallicity [dex] 


[Fe/H] 


+0.14 


±0.09 


M 


Luminosity [L©] 




4.25 


±0.41 


LC + M + T + Y^ 


Visual Magnitude [mag] 


Mv 


3.19 


±0.12 


LC + M + T + Y^ 


Apparent Magnitude [mag] 


my 


11.592 


±0.004 


my 


Age [Gyr] 




2.9 


±0.3 


LC + M + T + Y2 


Distance [pc] 




479 


±26 


LC + M + T + Y^ + my 


Planetary parameters: 










Mass [Mjup] 


M,, 


0.917 


±0.070 


LC + M + T+Y^ + K 


Radius [Rjup] 


Rp 


1.706 


±0.056 


LC + M + T + Y2 


Surface gravity [cgs] 




2.893 


±0.042 


LC + K 


Mean density [g cm"-'] 


pp 


0.229 


±0.027 


LC + M + T + Y2 



" The type of data used to calculate the given quantity. LC denote s light curve data. M and T denote metalUcity and 
temperature respectively and were taken from lSozzetti et all (l2009h . K, PRR, and my denote velocity semiampl itude, 
projected rotation rate, and apparent magnitude respectively and were taken from iMandushev et al.l (120071) . Y^ 
denotes the Yonsei-Yale isochrones (.Yi et aL 2001.) 



TABLE 6 
System Parameters for HAT-P-3 



Parameter 


Symbol 


Value 


68.3% Conf. Limits 


Comment" 


Transit parameters: 










Orbital period [d] 


P 


2.8997360 


±0.0000020 


LC 


Midtransit time [HJDujc] 




z45485o.701 18 


1 r\ r\c\{\ 1 c\ 

±0.00010 


LC 


Planet-to-star radius ratio 


J? / J? 

I 




±u.uu/u 




Orbital inclination [deg] 


/ .u / 


±U. JJ 




Scaled semimajor axis 


a/R-, 


10.39 


±0.49 


LC 


Transit impact parameter 


b 


0.530 


±0.075 


LC 


Transit duration [hr] 




2.075 


±0.022 


LC 


Transit ingress or egress duration [hr] 




0.270 


±0.033 


LC 


Other orbital parameters: 










Orbital eccentricity 


e 







Adopted 


Velocity semiamplitude [m s"' ] 


K 


89.1 


±2.0 


K 


Planet-to-star mass ratio 


Mp/M^ 


0.000615 


±0.000015 


LC + M + T+Y^ + K 


Semimajor axis [AU] 


a 


0.03866 


±0.00041 


LC + M + T + Y2 


Stellar parameters: 










Mass [Mq] 




0.917 


±0.030 


LC + M + T + Y2 


Radius [Rq] 


R* 


0.799 


±0.039 


LC + M + T + Y2 


Mean density [g cm"^] 


P* 


2.53 


±0.36 


LC 


Errective temperature [K] 


Teff 


5185 


±80 


T 


Projected rotation rate [km s"^] 


vsini* 


0.5 


±0.5 


PRR 


Surface gravity [cgs] 


loss* 


4.594 


±0.041 


LC + M + T + Y^ 


Metallicity [dex] 


[Fe/H] 


+0.27 


±0.08 


M 


Luminosity [L©] 




0.435 


±0.053 


LC + M + T + Y^ 


Visual Magnitude [mag] 


Mv 


5.87 


±0.15 


LC + M + T + Y^ 


Apparent Magnitude [mag] 


my 


11.577 


±0.067 


my 


Age [Gyr] 




1.6 


-1.3, +2.9 


LC + M + T + Y2 


Distance [pc] 




138 


±10 


LC + M + T + Y^ + my 


Planetary parameters: 










Mass [Mjup] 


M,, 


0.591 


±0.018 


LC + M + T+Y^ + K 


Radius [Rjup] 


Rp 


0.827 


±0.055 


LC + M + T + Y2 


Surface gravity [cgs] 




3.330 


±0.058 


LC + K 


Mean density [g cm"-'] 


Pp 


1.29 


±0.25 


LC + M + T + Y2 



" The type of data used to calculate the given quantity. LC denotes light curve data, K, M, T, and PRR de note velocity 
semiam plitude, metallicity, temperature, and projected rotation rate res pectively an d were taken from iTorres et al.l 
(l2007h . The uncertanties in M and T were adjusted as suggested by ^Torre s et al.l (120081) . my denotes apparent 
magnitude and was taken from .Droege et al.. (.2006,) . Y^ denotes the Yonsei-Yale isochrones (lYi et al...2001.) . 
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TABLE 7 
System Parameters for WASP- 12 



Parameter 


Symbol 


Value 


68.3% Conf. Limits 


Comment" 


Transit parameters: 










Orbital period [d] 


P 


1.0914222 


1 f\ r\r\f\f\f\ 1 1 

±0.0000011 


LC 


Midtransit time [HJDujc] 




2454500.97605 


±0.00028 


LC 


Planet-to-star radius ratio 


J? / J? 

I 


n 1 1 1 Q 


±u.uu/u 




Orbital inclination [deg] 


60. u 


±J.U 




Scaled semimajor axis 




3.097 


±0.082 


LC 


Transit impact parameter 


b 


0.22 


±0.15 


LC 


Transit duration [hr] 




3.001 


±0.037 


LC 


Transit ingress or egress duration [hr] 




0.324 


±0.020 


LC 


Other orbital parameters: 










Orbital eccentricity 


e 







Adopted 


Velocity semiamplitude [m s"'] 


K 


226 


±4 


K 


Planet-to-star mass ratio 


Mp/M^ 


0.000993 


±0.000038 


LCh-Mh-Th-Y^h-K 


Semimajor axis [AU] 


a 


0.02293 


±0.00078 


LC H- M H- T H- Y2 


Stellar parameters: 










Mass [Mq] 




1.35 


±0.14 


LC H- M H- T H- Y2 


Radius [Rq] 


R* 


1.599 


±0.071 


LC H- M H- T H- Y2 


Mean density [g cm"^] 


P* 


0.472 


±0.038 


LC 


Errective temperature [K] 




6300 


±150 


T 


Projected rotation rate [km s"'] 


vsiniV 


<2.2 


±1.5 


PRR 


Surface gravity [cgs] 


lege* 


4.162 


±0.029 


LC H- M H- T H- Y^ 


Metallicity [dex] 


[Fe/H] 


0.30 


±0.10 


M 


Luminosity [Lq] 




3.0 


±1.2 


LC H- M H- T H- Y^ 


Visual Magnitude [mag] 


Mv 


3.54 


±0.38 


LC H- M H- T H- Y^ 


Apparent Magnitude [mag] 


my 


11.69 


±0.08 


my 


Age [Gyr] 




1.7 


±0.8 


LC H- M H- T H- Y^ 


Distance [pc] 




427 


±90 


LC H- M H- T H- Y2 H- my 


Planetary parameters: 










Mass [Mjup] 


M,, 


1.404 


±0.099 


LCh-Mh-Th-Y^h-K 


Radius [Rjup] 


Rp 


1.736 


±0.092 


LC H- M H- T H- Y2 


Surface gravity [cgs] 




3.066 


±0.031 


LCh-K 


Mean density [g cm"-'] 


pp 


0.337 


±0.039 


LC H- M H- T H- Y2 



" The type of data used to calculate the given quantity. LC denotes light curve data. K, M, T, PRR, and my denote 
velocity semiam plitude, metallicity , temperature, projected rotation rate and apparent magnitude respectively, and 
were taken from iHebb et al.l (l2009h (the uncertainties were symmetrized). Y^ denotes the Yonsei-Yale isochrones 
JVi et al 2001.) . 



TABLE 8 
MiDTRANSiT Times for TrES-4 



Epoch Midtransit Time [HJD] Uncertainty Filter Source 



-2 


2454223.797215 


0.000847 


z 


Sozzetti et al. 


(2009) 





2454230.904913 


0.000656 


z 


Sozzetti et al. 


(2009) 





2454230.905624 


0.001106 


B 


Sozzetti et al. 


(2009) 


112 


2454628.947249 


0.001056 




This paper'' 




195 


2454923.920285 


0.000492 




This paper 




204 


2454955.907886 


0.001138 




This paper** 




305 


2455314.853547 


0.000616 




This paper 




307 


2455321.959132 


0.001049 




This paper** 





" Transits were o bserved by iMandushev et al.l (l2007h but re-analyzed by 
iSozzetti et all (12009') 
Partial transit observed 



TABLE 9 
Midtransit Times for HAT-P-3 



Epoch 


Midtransit Time [HJD] 


Uncertainty 


Filter 


Source 


-220 


2454218.759400 


0.002900 


i 


Torres et al. (2007) 


-111 


2454534.830232 


0.000391 


i 


This paper 


-101 


2454563.828281 


0.000358 


i 


This paper 


-91 


2454592.824589 


0.000707 


i 


This paper^ 





2454856.701370 


0.000240 


_b 


Gibson et aL (2010) 


1 


2454859.600240 


0.000370 


_b 


Gibson et aL (2010) 


17 


2454905.996437 


0.000385 


Z 


This paper 


19 


2454911.795918 


0.000510 


Z 


This paper^ 


22 


2454920.495670 


0.000550 


_b 


Gibson et aL (2010) 


28 


2454937.893861 


0.000317 


Z 


This paper 


31 


2454946.592600 


0.000650 


_b 


Gibson et al. (2010) 


32 


2454949.493340 


0.000400 


_b 


Gibson et al. (2010) 


41 


2454975.590370 


0.000340 


_b 


Gibson et al. (2010) 


42 


2454978.489930 


0.000510 


_b 


Gibson et al. (2010) 


Partial transit observed 








A single wideband filter (« 500 


-700nm) was used. 






TABLE 10 








Midtransit Times for WASP- 


12 


Epoch 


Midtransit Time [HJD] 


Uncertainty 


Filter 


Source 





2454508.97610 


0.00020 




Hebb et al. (2009)" 


304 


2454840.76781 


0.00047 


V 


This paper 


608 


2455172.56103 


0.00044 


V 


This paper 



" This data point actually represents a consensus value derived from observa- 
tions of several transits. 



